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For an unsteady baryonic gamma-ray burst (GRB) outflow, the fast and slow proton shells collide 
with each other and produce energetic soft gamma-ray emission. If the outflow has a significant 
neutron component, the ultra-relativistic neutrons initially expand freely until decaying at a larger 
radius. The late time proton shells ejected from the GRB central engine, after powering the regular 
internal shocks, will sweep these /3— decay products and give rise to very bright UV/optical emission. 
The naked-eye optical flash from GRB 080319B, an energetic explosion in the distant universe, can 
be well explained in this way. 

PACS numbers: 98.70.Rz 



There were four gamma-ray bursts (GRBs) detected 
on 19 March, 2008. Among them, GRB 080319B was 
most noticeable due to its huge isotropic energy and its 
extremely bright prompt optical emission that could have 
been seen with naked eye [l[. For a redshift z = 0.937, 
the detected peak optical emission had a visual mag- 
nitude 5.3 that corresponds to an optical luminosity 
> 5 X 10^° erg s~^. The simultaneous soft 7— ray emis- 
sion had a luminosity L-y ~ 4 x 10^^ erg s~^. The optical 
flux is significantly above the spectral extrapolation of 
the GRB emission into the optical band 1], suggesting 
that the optical and the gamma-ray emissions originate 
from different emission sites or belong to different spec- 
tral components at a same emission site. The lightcurves 
in the two bands show some similarities, but do not trace 
each other exactly [Ij. This suggests that the emissions 
from the two bands are somewhat related to each other. 

; For the energetic outflow of GRB 080319B, if not 
, strongly magnetized, a significant neutron component 
is unavoidable 0, S 0i [J- The average Lorentz fac- 
tor of the outflow before getting decelerated by a stel- 
lar wind medium is very high. A lower limit can be 
set by the Lorentz factor of the forward shock at ~ 70 
s, when the X-ray afterglow began to decline normally, 
i.e., r > 500i;i[^t;^;,'_4^V70s)-i/4[(l + z)/2]i/4, where 
i?k ^ 10^^ erg is the isotropic-equivalent energy of the 
outflow and ^ 0.01 is the stellar wind parameter of 
the progenitor [ij. The current data then suggests an 
initial Lorentz factor of the outflow F ~ 10'^. For an 
unsteady baryonic outflow, the GRB is powered by the 
interaction of proton shells with variable Lorentz factors, 
i.e., internal shocks [1, 0, H, 0, • In order to convert a 



significant fraction of the initial kinetic energy into inter- 
nal energy and then to 7— ray radiation, the difference in 
Lorentz factors between the shells should be substantial 
(i.e., the Lorentz factors ?7f 3> ?7s, hereafter the subscript 
"f " and "s" denote "fast" and "slow" , respectively) and 
their masses should satisfy Aff = /Mg with an / > rjs/rji. 
The mergered new proton shell moves with a Lorentz fac- 
tor F„i « ^/Jrjfrk S]. For ?7f ~ a few x iC and / ~ 0.1, 
^ a few X 100 is needed to get a rm(~ F) ^ 10^. 

The dynamics of a neutron-rich outflow is governed by 
the dimensionless entropy rj = L/ [Mc^) at tq, where M is 
the mass loading rate, L is the isotropic luminosity of the 
ejecta, and rg is the radius of the central engine. Whether 
or not the proton and the neutron components decouple 
from each other depends on whether 77 is above or be- 
low the critical value rycr ^ lO^Ll{^r'^lj'\\ + Cl/'A'^''^, 
where ^ is the ratio of the number density of neutrons 
to protons The convention Qx = Q/10^ is adopted 



in this work in cgs units. For rj < rjcr, the neutron and 
proton components (denoted by the subscripts n and p, 
respectively) are still coupled with each other by nuclear 
elastic scattering at a radius ~ Tyro, i.e., at the end of the 
flreball acceleration. Thus for slow shells one usually has 
??n,s = Vp,s- For 77 > ?7cr, the n and p components decou- 
ple at a radius ~ rjnjro where the nuclear elastic scatter- 
ing becomes too weak to accelerate the neutrons, where 

For the observed information of GRB 080319B, we take 
ilp,s — Vn,s{^ Vnj) ^ a few hundred, rjpj ~ a few x 10^, 
Afp.s ~ lOMp^f, and ^ ^ 1. Hereafter Fn denotes the 
generic Lorentz factor of the neutron shell regardless of 
whether it is relatively "slow" or "fast" . Below we show 
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TABLE I: Physical parameters that reproduce the prompt 
emission of GRB 080319B. 



Quantity 


slow shell 


fast shell 


Lorentz factor of protons 


~ 400 


~ 2 X 10" 


Lorentz factor of neutrons 


~ 400 


~ 500 


isotropic luminosity (erg/s) 


- 2 X lO-'^^ 


~ 10^* 


mass ratio of neutrons to protons 


~ 1 


~ 1 


entropy per baryon {s/k) 


~ 1.2 X 10^ 


~ 4 X lO** 



that with these parameters, both the soft 7— ray and op- 
tical emission can be interpreted. 

Firstly we discuss the regular internal shocks powered 
by the collisions of the fast and slow proton shells at a ra- 
dius i?i„t ~ 2Vl^c5t/{l+z) - 3x10" cmr^ 3(5t_2[2/(l-h 
z)] cm, where 5t is the typical variability timescale of the 
prompt 7— ray light curve. The electrons are accelerated 
by the internal reverse shock to a typical Lorentz factor 
d ~ 5 X 103(ee/0.3)(r7pj/2 x 10'')r-|3 and the mag- 
netic field generated in the shocks can be estimated as 
^ 4x 103 Gauss {2,eB/e,Y'^L\% ,,R7-^,^^^ ,;C^\, where 
and es are the fractions of the shock energy distributed 
to electrons and magnetic fields, respectively. The typi- 
cal synchrotron emission frequency of these electrons in 
the internal reverse shock is ~ 1.4 x 10^° Hz, matching 
the observation. The internal forward shock can only ac- 
celerate electrons to a typical Lorentz factor ^ 300 and 
the emission is in the soft X— ray band. 

The neutrons have negligible interaction with the pro- 
tons before decaying into protons, electrons and electron 
neutrinos. The /3— decay radius reads 



Re 



1.1 X 10^*^ cm r 



n,2.6- 



(1) 



A pair of fast/slow proton shells ejected at late times 
would merge into a proton shell with a Lorentz factor 
^ Fin in the inner internal shocks, and then catch up 
with the decay trail of the neutron shells ejected ear- 
lier at a radius i?cat ~ 2F^C(5T/(1 -I- z), where 5T is the 
ejection time-lag between the earlier neutron shell and 
the later proton shell. As long as i?cat > ^/3, which re- 
quires 5T > 1.1(1 -I- 2)F~2.6 sec, there will be a substan- 
tial amount of /?— decay products that would be swept 
orderly by the later proton shell at a radius of ~ 2R/s. 
These "secondary internal shocks" are unable to give rise 
to energetic X— ray and 7— ray emissions for the following 
reasons. (1) They are generated at a radius much larger 
than i?int, so that the magnetic fields are much weaker 
than those in the regular internal shocks. (2) They have 
much lower efficiency than the inner ones at i?int because 
of the smaller Lorentz contrast between the merged pro- 
ton shell and the neutron shell. For the same reason, the 
total energy converted into internal energy and then into 
radiation is lowered by a factor of > 10 than that of the 
regular internal shocks. 

We first look at the interaction between a proton shell 
formed in the inner internal shocks and the decay trail 



of a series of identical neutron shells. The number den- 
sity of the decay products being swept by the proton 
shell reads [131 n ~ p2^"p , where m,, is the neu- 
tron rest mass, and M„ oc exp(— is the rest mass 
of the neutron shell that undergos decay. The mini- 
mum Lorentz factor of the shocked electrons can be es- 
timated as y 7e,m ~ 55 ( ^ ) ^IpZl] ( ^^"f^ ) , where p is 
the power-law index of the accelerated electrons, and 
7rci ~ (r/rn + r,i/F)/2 is the Lorentz factor of the trail- 
ing fast-moving proton shell with a Lorentz factor F rel- 
ative to the neutron shell. The magnetic field strength 
B' can be estimated as fl^ 

B' ^ 60 Gauss ee^-iCTrci " 

Ri6R^\6^n,L exp{-R/2Rp), (2) 

where the neutron number Nn of one shell is estimated by 
Nn ^ ^L6t/[{1 + z)T^mnC^] ~ 3.3 x lO^i ^L5iSt-2r-%. 
Since this secondary internal shock region is permeated 
by gamma-ray photons produced from the "inner" inter- 
nal shocks from the late-time ejected proton shells, the 
electrons in the secondary internal shock reg ion suffer 
Compton cooling by these prompt 7— rays [1411. The cor- 
responding cooling Lorentz factor reads 



7e,c 180 Fgi^igi^ 52.7- 



(3) 



So, the synchrotron radiation energy of the shocked 
electrons peaks at a frequency 

/.p ^ 1.4 X IQi^ Hz r^jg. "^M%%^7e^,n.} _ 

The detected maximum specific spectral fiux can be 
estimated as llSl 



F^n,.. - 50 Jy A^e.ssraSz, 



(5) 



where - NJtp/St ~ 2 x lO'^^ ^L^iT-^^^T-^Rp^ie, and 
Stf3 « (1 z)Ri3/r'^c ~ 0.6 s (i±^)i?^,i6F-2. This flux is 
bright enough to well exceed the spectral extrapolation 
of the gamma-ray emission and to interpret the observed 
naked-eye optical flash. Notice that we have introduced 
the number of proton shells (~ Stp/St) ejected during 
the time span of Stp to account for the total emission 
output tia]. Since F^^^^^ cx NeB' cx ^^/^^ ^ ~ 1 is highly 
needed to reproduce the prompt optical emission with a 
flux - 20 Jy. 

Following the standard approach [l^ , we estimate the 
synchrotron-self absorption frequency 



2 X 10" Hz 



100 



min{7e,c,7e,Tn} 



,.3/5 p-6/5 p/2/5 
^^n, sr. 5-^^16 



s';-F3. 



(6) 

We can see that for the standard parameters, in partic- 
ular the large neutron decay radius Riq ^ 1 and the 
small -B2 < 1, one has < i^p- This is the main rea- 
son that the optical flash of this burst is so bright. The 
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lack of bright optical flash in most other bursts may be 
attributed to their smaller Lorentz factors, which give a 
smaller neutron decay radius, stronger magnetic fields, 
and hence, a higher self-absorption frequency than opti- 
cal. In addition, smaller F and Rp would suppress the 
optical emission also by reducing 7e,c, as shown in eq.([n]). 

Till now we have shown that with reasonable parame- 
ters (see TabUfor a summary), the neutron-rich internal 
shocks can power both the energetic soft 7— ray flare and 
the extremely bright optical flash of GRB 080319B. The 
thermodynamic entropy per baryon of the initial ejecta 
can be estimated as s//c ~ rjm^c^ /{kT) [3,[lH|) where k is 
the Boltzman's constant and T ~ 5 MeV r^y^Ly^ is the 
temperature of the initial flreball. For the fast and slow 
shells, we have s/k ~ (4 x 10^, 1.2 x 10^), respectively. 

Since the neutron shells are originally coupled to the 
early proton shells, they carry the essential variability 
information of the early proton outflows. When the de- 
cay products are swept by a later injected proton shell, 
the resultant optical lightcurve generally follow the vari- 
ability pattern of the gamma-ray lightcurve, yet does not 
strictly trace the gamma-ray lightcurve. The time delay 
between the optical and gamma-ray peaks can be esti- 
mated as ~ 2.2[(l-|-z)/2]F~2 g ^^'^i matching the time lag 

3 sec found in the correlation analysis of prompt 7— ray 
and optical lightcurves Being at a larger emission 

radius ~ Rp^ the optical variability is also smoothed by 
the geometric effect in a timescale ~ 5tp ^ 5t, gener- 
ally consistent with the fact that the optical lightcurve 
is much smoother than the gamma-ray lightcurve [3]. 
The contamination from a bright external reverse shock 
optical emission component jH introduces a bright back- 
ground and also contributes to the smoothing. 

In the s ync hrotron self-Compton (SSC) internal shock 
model p. ]l7t. one requires a large emission radius {'^ 
10^^ cm) to avoid self-absorption, and also requires a 
large Ce/es (> 10^) to interpret the huge energy contrast 
between the optical and the 7— ray emission. In that sce- 
nario, one predicts a very bright prompt GeV emission 
component with luminosity > 4 x 10^'^ erg s~^, which 
may suffer an energy crisis jl9| . In our scenario, the SSC 
component of the prompt GeV emission is most likely 
weaker because our standard model parameters demand 
te cb. GeV emission can also arise from the prompt 
7— ray cooling in the shocked region at ~ i?^. However, 
as mentioned before, the total energy of the shocks at 
~ Rf3 is much smaller than that of the "inner" inter- 
nal shocks. So the high energy emission would not be 
enhanced significantly. Prompt UV/optical photons will 
cool the forward shock electrons effectively and give rise 
to GeV-TeV emission with Luminosity ~ 10^^ erg s^^. 
However, such a component, appearing as a plateau, 
should last two more times longer than the prompt opti- 
cal emission and thus can be easily distinguished 20]. 

In summary, we show in this work that an outflow, 
containing comparable amounts of protons and neutrons 



(i.e., '^'^ 1) and having an averaged Lorentz factor ~ 10^ 
and a huge isotropic luminosity L ~ 10^^ erg s~^, can 
power an energetic soft 7— ray burst and a naked-eye op- 
tical flash in GRB 080319B (see Tab|l!]for a summary). 
In our scenario, the prompt optical light curve will be 
much smoother than that of the prompt 7— ray emission 
and the peaks in optical band will lag behind the 7— ray 
peaks by a few seconds, both are well consistent with 
the data 17, [3] and are hard to be interpreted in 



other models that invoke a similar emission radius for 
both prompt gamma-ray and optical emissions P.[l7. 21|. 
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TABLE II: General features of the neutron-rich internal shock model for GRB 0803 19B. 



Quantity 



regular internal shocks 



secondary internal shocks 



slow material 
shock radius 
shock strength 
typical emission 
emission flux 



slow protons from central engine 
~ 10" - 10^^ cm 
ultra-relativistic 
7— rays 

~ lO"'' erg/s/cm^ 



variability timescale 3t ~ 0.01 — 0.1 sec 



/3— decay products of neutrons 

Rf) ^ 10^^ cm 
sub-relativistic 
ultraviolet/optical photons 
~ 50 Jy (in optical band) 



1 sec 
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